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Nonthermal radiation observed from astrophysical systems containing relativistic jets and shocks, e.g., active galactic 
nuclei (AGNs), gamma-ray bursts (GRBs), and Galactic microquasar systems usually have power-law emission spectra. 
Recent PIC simulations using injected relativistic electron-ion (electro-positron) jets show that acceleration occurs 
within the downstream jet. Shock acceleration is a ubiquitous phenomenon in astrophysical plasmas. Plasma waves and 
their associated instabilities (e.g., the Buneman instability, other two-streaming instability, and the Weibel instability) 
created in the shocks are responsible for particle (electron, positron, and ion) acceleration. The simulation results show 
that the Weibel instability is responsible for generating and amplifying highly nonuniform, small-scale magnetic fields. 
These magnetic fields contribute to the electron's transverse deflection behind the jet head. The "jitter" radiation 
from deflected electrons has different properties than synchrotron radiation which assumes a uniform magnetic field. 
This jitter radiation may be important to understanding the complex time evolution and/or spectral structure in 
gamma-ray bursts, relativistic jets, and supernova remnants. 
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1. Introduction 

The non-thermal spectrum of GRB sources is 
thought by most to arise in shocks which develop 
beyond the radius at which the relativistic fireball 
has become optically thin to 77 collisions [1]. How- 
ever, the observed spectra are hard, with a signifi- 
cant fraction of the energy above the 77 — ► e + e~ 
formation energy threshold. A high compactness pa- 
rameter can result in new pairs being formed out- 
side the originally optically thin shocks responsible 
for the primary radiation [2] . Radiation scattered by 
the external medium, as the collimated 7-ray front 
propagates through the ambient medium, would be 
decollimated, and, as long as I > 1 absorbed by the 
primary beam. An pair cascade can then be pro- 
duced as photons are back-scattered by the newly 
formed pairs and interact with other incoming 
seed photons [2-8]. Here, I = Lctt / Aitr\m, a (? , where 
L is the total luminosity, and r r is the characteris- 
tic source dimension. In this paper, we consider the 
plasma instabilities generated by rapid e^ 1 pair cre- 
ation in GRBs. The injection of pairs induces 
strong streaming motions in the ambient medium 



ahead of the forward shock. This sheared flow will 
be Weibel-like unstable [9], and if there is time be- 
fore the shock hits, the resulting plasma instabili- 
ties will generate sub-equipartition quasi-static long- 
lived magnetic fields on the collisionless temporal 
and spatial scales across the pair-enriched region. 
This process is studied using three-dimensional ki- 
netic simulations of monoenergetic and broadband 
pair plasma shells interpenetrating an unmagnetized 
medium [10]. 



2. Parameters and Initial Conditions 
used in Recent Simulations 

Simulations have been performed using an 85 x 85 x 
640 grid with a total of 380 million particles (27 
particles/cell/species for the ambient plasma) and an 
electron skin depth, A cc = c/w pc = 9.6A, sufficient 
to study nonlinear spatial development [11-13]. The 
time step t = 0.013/wpe where w pc = (e 2 n c /rn c ) 1 ^ 2 
is the electron plasma frequency (n e = rib where rib 
is the ambient "background" plasma density). 
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Fig. 1. Snapshots viewed from the front of the jet at t 



59.1 



left panel: the isosurface of Z-component of the current 



density (±J Z ) with the magnetic field lines (white) at the linear stage, 
current density (— J z ) with the magnetic field lines behind the jet front. 



and right panel: the isosurface of Z-component of the 



3. Monoenergetic &: Broad Pair Jet 
distributions Injected into Pair & 
Electron-Ion Plasmas 

The simulations described below were performed by 
a newly parallelized OpenMP code on Columbia at 
NASA Advanced Supercomputing (NAS). 

We have simulated four different initial pair 
jet distributions. Case A is a mono-energetic jet 
(7V1I = 12.57c) injected into a pair ambient plasma 
and similar to published simulations [11, 12]. Case 
B is also a mono-energetic jet but injected into an 
electron-ion ambient plasma. Cases C and D corre- 
spond to pair jets with a broad velocity distribu- 
tion (2 < Mm = 'yVjj/c < 30) created by photon 
annihilation [5]. In case C the jet is injected into 
a pair ambient plasma and in case D the jet is in- 
jected into an electron-ion ambient plasma. For all 
cases the jet particles are very cold (0.01c in the 
rest frame). The average jet velocity for the broad 
pair loaded distribution is 'jVu/c ~ 12.5. The mass 
ratio of electrons to ions in the ambient plasma is 
mi/m e = 20. The electron thermal velocity in the 
ambient plasma is i»th = 0.1c, where c is the speed of 
light. The ion thermal velocity in the ambient plasma 
is Vthi — 0.022c. Otherwise parameters and setup are 
the same as those in the previous simulations. 

Encountering the medium at rest, incoming 
pairs are rapidly deflected by field fluctuations. The 
initial perturbations grow nonlineary as the deflected 
pairs collect into current channels. The resultant 
toroidal magnetic fields cause mutual attraction be- 
tween currents forcing like currents to approach each 
other and merge. As a result, the magnetic field 



grows in strength. This continues until the fields grow 
strong enough to deflect the much heavier ions [14]. 
The ions stay clearly separated in phase space and 
are only slowly heated. In the presence of ions, the 
magnetic field saturates at a higher level by a fac- 
tor of (mi/m,,) 1 / 2 = \/20 ~ 4.5, albeit on a longer 
timescale. In electron-ion jets, the real mass ratio 
of 1836 will be expected to provide an even higher 
level of magnetic field and enhanced electron accel- 
eration [15, 16]. 

The differences between an ambient consisting of 
pairs and an ambient consisting of electrons and ions 
are due to the massive ion bulk momentum constitut- 
ing a vast energy reservoir for particle heating [10]. 
We also found that a broad pair jet velocity distribu- 
tion sustains a stronger magnetic field over a longer 
region than a comparable mono-energetic pair distri- 
bution. This additional effect is a result of the more 
"massive" high Lorentz factor pairs in the jet. 

The electron density and current filaments have 
a complicated three-dimensional structure. Current 
filaments (J z ) and associated magnetic fields (white 
curves) produced by the filamentation (Weibel) in- 
stability form the dominant structures in the rela- 
tivistic collisionless shock shown in Figure 1. In the 
linear stage, the transverse size of these structures 
is nearly equal to the electron skin depth but the 
longitudinal size (along the jet direction as shown in 
the left side of the left panel) is much longer. The 
growing current filaments in the linear stage merge 
in the nonlinear stage (behind the jet front). This 
appears spatially as merging behind the jet front of 
the smaller scale filaments that appear first. 
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Fig. 2. The path of a charged particle moving in a homogenous magnetic field (left panel). The particle radiates a time dependent 
electric field. An observer situated at great distance along the n-vector sees the retarded electric field from the gyrating particle 
(right panel). As a result of relativistic beaming, the field is seen as pulses peaking when the particle moves directly towards the 
observer. 



4. Radiation from a gyrating electron 

Let a particle be at position r (t) at time t (see 
Fig. 5 in Hededal [15]). At the same time, we ob- 
serve the electric field from the particle from posi- 
tion r. However, because of the finite propagation 
velocity of light, we observe the particle at an ear- 
lier position ro (t ) where it was at the retarded time 
t = t-St = t- R(t')/c. Here R(t) = |r - r (t')| 
is the distance from the charge (at the retarded time 
t ) to the observer point. 

The retarded electric field from a charged parti- 
cle moving with instant velocity (3 under acceleration 
$ is expressed as Jackson [17]. 
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Here, n = R(t )/|R(t )| is a unit vector that 
points from the particles retarded position towards 
the observer. The first term on the right hand side, 
containing the velocity field, is the Coulomb field 
from a charge moving without influence from exter- 
nal forces. The second term is a correction term that 
arises, if the charge is subject to acceleration. Since 
the velocity-dependent field is falling off as R~ 2 while 
the acceleration-dependent field falls off as R -1 , the 
latter becomes dominant when observing the charge 
at large distances (R ^> 1). 

After some calculation and assumptions (for a 
detailed derivation see e.g. Hededal [15]) the total 
energy W radiated per unit solid angle per unit fre- 



quency is then expressed as 
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This equation contains the retarded electric field 
from a charged particle moving with instant veloc- 
ity (3 under acceleration (3, and only the accelera- 
tion field is kept since the velocity field decreases by 
1/R 2 . The distribution over frequencies of the emit- 
ted radiation depends on the particle energy, radius 
of curvature, and acceleration. These quantities are 
readily obtained from the trajectory history of each 
charged particle. 




Fig. 3. The observed power spectrum from a single charged 
particle, gyrating in a magnetic field with different viewing 
angles. The viewing angles are 0°, 1°, 2°, 3°, 4°, 5°, and 6° 
(n y 7^ 0) and their peak frequencies are 7642, 4395, 1648, 666, 
316, 166, and 133, respectively. With larger angles the fre- 
quencies above the Nyquist frequency arc strongly distorted. 
The units on the axis arc arbitrary. 

The radiation from single electron gyrating in 
the uniform magnetic field (Bq = B y ) is calculated 
as shown in Fig. 2. The Lorenz factor of the electron 
is 40.8. The spectra observed far from the electron 
in the z direction are shown in Fig. 3. The critical 
frequency f c = | 7 3 (j^ = 2309, where p = 11.03. 
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The higher frequencies (> / c ) are strongly damped 
with increasing angles as in Jackson [17]. 

5. Discussion 

We present self-consistent three-dimensional simula- 
tions of the fields developed by the electromagnetic 
filamentation instability as a consequence of e^ 1 pair 
injection, which arises naturally in the environment 
of GRB sources as a result of back-scattering. Our 
results demonstrate that even in an initially unmag- 
netized scattering plasma, a small-scale, fluctuating, 
predominantly transversal, and near-equipartition 
magnetic field is unavoidably generated. These fields 
maintain a strong saturated level on spacial scales 
much longer than A c at least for the duration of 
the simulations ~ 60(c/A o ) <~ 0.1(c/At) • The e 1 * 1 
pairs are effectively scattered with the magnetic field, 
thus effectively communicated their momentum to 
the scattering medium initially at rest. Our results 
indicate that the fields necessary to ensure that e^ 1 
pairs remain coupled to the medium can be easily 
created via plasma instabilities. The next required 
step is to increase the ion-to-electron mass ratio in 
the scattering medium in order to determine the spa- 
tial spread and character of the particle coupling 

In order to obtain the spectrum of synchrotron 
(jitter) emission in future work [18-21], we will con- 
sider an ensemble of selected electrons in the region 
where the Weibel instability has grown fully and 
electrons are accelerated in the generated magnetic 
fields. We will calculate emission from about 50,000 
electrons during the sampling time, t s = ti — t\ 
with Nyquist frequency wn = 1/2 At where At is 
the simulation time step and the frequency resolu- 
tion Alu = l/t s . However, since the emission coordi- 
nate frame for each particle is different, we accumu- 
late radiation at a fixed angle in simulation system 
coordinates after transforming from the individual 
particle emission coordinate frame. This provides an 
intensity spectrum as a function of angle relative to 
the simulation frame Z-axis [can be any angle by 
changing the unit vector n in eq. (1)]. A hypothet- 
ical observer in the ambient medium (viewing the 
external GRB shock) views emission along the sys- 
tem Z-axis. This computation is carried out in the 
reference frame of the ambient medium. For an ob- 
server located outside the direction of bulk motion, 
e.g., internal jet shocks or an external relativistic jet 
driven bow shock, an additional Lorcntz transforma- 



tion is needed along the line of sight to the observer. 
The spectra obtained from the simulations need to 
be rescaled with a realistic time scale and relativistic 
Doppler shift. In the electron-ion jets, the real mass 
ratio 1836 will provide enhanced electron accelera- 
tion compared to a mass ratio of 20 [15, 16]. 
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